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Abstract We propose to characterize heavy-ion colli¬ 
sions at ultra-relativistic energies by using fluctuations 
of energy density and temperature. Temperature fluc¬ 
tuations on an event-by-event basis have been stud¬ 
ied both in terms of global temperature of the event, 
and locally by constructing fluctuation maps in small 
phase space bins in each event. Global temperature fluc¬ 
tuations provide an estimation of the specific heat of 
the system. Local temperature fluctuations of the event 
may be ascribed to the remnants of initial energy den¬ 
sity fluctuations. Together these two observables give 
an insight into the system created in heavy-ion colli¬ 
sions and its evolution. Event-by-event hydrodynamic 
calculations indeed provide adequate theoretical basis 
for understanding the origin of the fluctuations. We 
demonstrate the feasibility of studying global and local 
temperature fluctuations at the Large Hadron Collider 
energy by the use of AMPT event generator. 
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1 Introduction 

Heavy-ion collisions at ultra-relativistic energies cre¬ 
ate matter at extreme conditions of energy density (e) 
and temperature (T), similar to the ones that existed 
within a few microseconds after the Big Bang [I]. The 
fireball produced in the collision goes through a rapid 
evolution from an early partonic phase of deconfined 
quark-gluon plasma (QGP) to a hadronic phase and ul¬ 
timately freezing out after a few tens of fm. The major 
goals of colliding heavy-ions at the Relativistic Heavy 
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Ion Collider (RHIC) and at the Large Hadron Collider 
(LHC) are to study the nature of the phase transition 
and understand the QGP matter in detail. With the 
production of large number of particles in each of the 
collisions, it has become possible to extract thermo¬ 
dynamic quantities on an event-by-event basis, rather 
than averaging over a sample of events. Event-by-event 
fluctuations of e, T, mean transverse momentum, par¬ 
ticle multiplicity, particle ratios, etc., as well as fluc¬ 
tuations of conserved quantities within finite detector 
acceptances have been proposed to provide dynamical 
information regarding the evolving system [2l[3l|4l[5j|6l 



In this manuscript, temperature fluctuations have 
been studied on an event-by-event basis by estimating 
global temperature of the event, and local temperatures 
in small phase space bins within the event. Determina¬ 
tion of temperature and its fluctuation for heavy-ion 
collisions have been possible because of large number 
of particles emitted in each event min], which is es¬ 
sential to keep the statistical uncertainties of the mea¬ 
surements within a reasonable limit. Global event to 
event temperature fluctuations provide the heat capac¬ 
ity (C), which is an important thermodynamic quantity 
in characterizing the system. Lattice QCD calculations 
predict a strong temperature dependence of the heat ca¬ 
pacity, the nature of which depends on the order of the 
phase transition ms\- This study sheds light on the 
nature of phase transition at the LHC energy. Fluctu¬ 
ations of initial energy density and temperature may 
survive till the freeze-out and manifest themselves in 
the measured temperature fluctuations. In fact, initial 
fluctuating conditions have been found to be necessary 
for explaining observed elliptic flow in central collisions 
and substantial triangular flow of charged particles m- 
The initial state fluctuations may have their imprint on 
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Fig. 1 (Color online). Distributions of energy density (upper panles) and temperature (lower panels) in the transverse (X-Y) 
plane at four proper times (t) obtained from hydrodynamic calculations for a single central Pb-Pb event at =2.76 TeV. 


the bin to bin local fluctuations within an event. Event- 
by-event hydrodynamic calculations provide a strong 
theoretical basis for studying the global and local tem¬ 
perature fluctuations. 


2 Event-by-Event hydrodynamical approach 


Local fluctuations in energy density arise because of 
the internal structures of the colliding nuclei. These ini¬ 
tial fluctuations manifest into local temperature fluctu¬ 
ations of the fireball at different stages of the collision. 
Relativistic hydrodynamic calculations which take such 
effects into account reveal large local fluctuations in e 
and T in small phase space bins at early stages of the 
collision [DIIIIITIIIH]. The local fluctuations have been 
quantified throughout the evolution by simulating cen¬ 
tral (0-5% of the total cross section) Pb-Pb events at 
^snn = 2.76 TeV by the use of a (2-|-l)-dimensional 
event-by-event ideal hydrodynamical framework m with 
lattice-based equation of state [19] . The formation time 
of the plasma is taken to be 0.14 fm irniiii]. A wounded 
nucleon (WN) profile is considered where the initial en¬ 
tropy density is distributed using a 2-dimensional Gaus¬ 
sian distribution function, 


s{X,Y) 


K 

27r(7^ 


NwN 

E 

2=1 


exp I 


(A - + (V 

2a2 



Here Xi , V are the transverse coordinates of the nu¬ 
cleon and K is an overall normalization constant. The 
size of the density fluctuations is determined by the free 


parameter cr, which is taken to be 0.4 fm |16] . The tran¬ 
sition temperature from the QGP to the hadronic phase 
is chosen to be 170 MeV and the kinetic freeze-out tem¬ 
perature is taken as 160 MeV. The results for e and T 
at different times (t) are analyzed for each collision in 
X-Y phase space bins in the transverse plane (each bin 
is chosen to be 0.6 fmxO.6 fm). 

Time evolutions of the distributions of e and T have 
been presented for four values of r in Fig. [T] for a sin¬ 
gle event. The upper panels (a-d) of the figure show 
three dimensional view of e, whereas the lower panels 
(e-h) show corresponding values of T in the transverse 
plane. At early times, sharp and pronounced peaks in e 
and hotspots in T are observed. Large bin-to-bin fluc¬ 
tuations observed in e and T indicate that the sys¬ 
tem formed immediately after collision is quite inho¬ 
mogeneous in phase space. As time elapses, the system 
cools, expands, and the bin-to-bin variations in e and 
T smoothen out. 

Observations from Fig. [T] have been quantified in 
terms of the mean energy density ((e)), mean tempera¬ 
ture ((T)) over the X-Y bins, and the bin-to-bin fluctu¬ 
ations of e and T as a function of r. The time evolution 
of (e), (T), and their fluctuations are presented in Fig.j^J 
where r is plotted in logarithmic scale for zooming in 
on the early times. The shaded regions represent the 
extent of event-by-event variations, taken from a sam¬ 
ple of five hundred events. It is observed that (e) falls 
sharply from ~168 GeV/fm^ at r = 0.14 fm to a value 
of ~20 GeV/fm^ at r = 1 fm, and then falls slowly 
till freeze-out. The initial e is close to the experimental 
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Fig. 2 (Color online). Temporal evolution of (a) average en¬ 
ergy density, (b) average temperature, (c) fluctuations in en¬ 
ergy density, and (d) fluctuations in temperature for central 
Pb-Pb collisions at y'SNN =2.76 TeV, obtained from hydro- 
dynamic calculations. Averaging is taken over the X-Y bins 
in every event. The shaded regions represent the extent of 
event-by-event variations for a large number of events. 


result, estimated by the ALICE collaboration [T31I21]- 
The fall of (T) with r is smooth, which goes down from 
^530 MeV at r = 0.14 fm to ~300 MeV at r = 1 fm. 
At the freeze-out, as expected, (T) is close to 160 MeV. 

The bin-to-bin fluctuations, Ae/(e) and AT/{T), 
are presented in the right panels of Fig. [5J where Ae 
and AT are the root mean square (RMS) deviations 
of the two quantities over the bins. At early times, the 
fluctuations are observed to be very large (^90% and 
^35% for e and T, respectively), which indicate the 
violent nature of the matter created in the collision. In¬ 
terestingly, although (e) decreases quite fast, the fluc¬ 
tuation in e remains almost constant up to r^2.5 fm, 
and then decreases rapidly. Around the same r, the 
fluctuation in T shows a kink, beyond which the fluctu¬ 
ation decreases even faster. During the hydrodynamic 
evolution, there may be a characteristic change in the 
behaviour of the system at this time, which needs to 
be understood. Nevertheless, it is clear that a detailed 
insight to the evolution of fluctuations is possible by 
studying local temperature fluctuations. 


3 Global Temperature Fluctuations 

Fluctuations of any observable of a system have two dis¬ 
tinct origins, first, quantum fluctuations which are ini¬ 
tial state fluctuations occurring at fast time scales and 
second, classical thermodynamical fluctuations which 
occur after elapse of sufficient time after the collision [531 


IMj- Initial state fluctuations arise because of internal 
structures of the colliding nuclei, and fluctuation in ini¬ 
tial energy densities, and appear as event-by-event fluc¬ 
tuations in the energy density or temperature. Ther¬ 
modynamic fluctuations have multiple sources, such as 
local thermal fluctuations of energy density, and event- 
by-event variation in the freeze-out conditions. Fluctua¬ 
tions of these observables are related to several thermo¬ 
dynamic parameters. For example, fluctuation in tem¬ 
perature are related to the heat capacity of the sys¬ 
tem wm- 


1 (AT^) 

c ~ r2 ’ 


( 2 ) 


where T is the temperature of the system at freeze-ou10, 
and Z\T2 = (T2) - (T)2. 

Fluctuations in a thermodynamic quantity like tem¬ 
perature may arise due to fluctuations in energy or en¬ 
tropy or finite particle number, and can be calculated 
from (pt) (mean transverse momentum) fluctuations 
or pt correlations However, ac¬ 

cording to Ref. |5], variation of C is similar to fluc¬ 
tuation in energy rather than the px fluctuation. Thus 
px correlation is not an appropriate measure of the tem¬ 
perature fluctuation, and may also affect measurement 
of C. To circumvent this limitation, we have come up 
with a new approach by assigning a temperature to an 
event and estimating event-by-event fluctuation for a 
given class of events. In fact, for central Pb-Pb collisions 
at the LHC energy, large number of produced particles 
in each event form a well defined exponential distribu¬ 
tion in Px. The inverse slope parameter of the distribu¬ 
tion gives the effective temperature, and subsequently, 
the heat capacity of the system can be calculated by 
using eqn. ©■ 

A feasibility study of calculating the global event 
temperature and local temperatures in small phase space 
bins has been made for Pb-Pb collisions at .^/snn = 2.76 
TeV, corresponding to the LHC energy. The string melt¬ 
ing (SM) mode of the AMPT model [35] has been em¬ 
ployed, which mimics the experimental conditions. This 
mode of AMPT model, using parameter set-B as given 
in Ref. [361I37] . includes a fully partonic phase that 
hadronizes through quark coalescence, and has been 
shown to reproduce the experimental data at LHC en¬ 
ergy [SSjlMlISZlISE]- 

A large number of particles produced in every cen¬ 
tral Pb-Pb collisions at ^snn =2.76 TeV makes it pos¬ 
sible to construct mx spectra of identified particles in 


^ For part of a thermodynamically equilibrated system with 
fluctuating energy {E), the heat capacity [^ can be expressed 
through (AE^) = T^C{T). This can be applied for a locally 
thermalized system produced during the evolution of heavy- 
ion collisions. 
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Fig. 3 (Color online). Event-by-event distribution of 
(mx) and corresponding Teff for pions within one unit of 
rapidity and full azimuth for central Pb-Pb collisions at 
y^snn = 2.76 TeV using AMPT. Distributions for synthetic 
events representing the statistical component are overlaid on 
the figure. 
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every event. Fitting the mx distribution of pions with 
a Maxwell-Boltzmann distribution yields the value of 
the effective temperature, Teg in each event. This tem¬ 
perature is related to the (mx) of pions within a given 
acceptance: 


(tot) 


2T^g + 2moTes + ml 

Too + Teff 


(3) 


where toq is the mass of pions. From the AMPT events, 
it is verified that for pions within the central rapid¬ 
ity (—0.5 < y < 0.5), full azimuth and toq < tox < 
1.5 GeV, the value of Teff from Htting the tox spectrum 
as well as from (tox) are similar to each other. The up¬ 
per limit of Tox is chosen to exclude pions, which may 
be affected by jets. Figure [3] shows the distributions 
of (tox) of pions and corresponding values of Teff cal¬ 
culated using the above expression for a large num¬ 
ber of events. It should be noted that the values of 
Teff have contributions from both a thermal part (Tkin) 
and a second component which depends on the collec¬ 
tive transverse velocity ((/3x)) of the system such that 
Teff = Tkin + /(/3t)- For pions, /(/3x) « mo{(3T)'^ [^. 
For top central (0-5%) collisions, the transverse velocity 
can be assumed to be same for all events due to spher¬ 
ical symmetry, and the fluctuations in the flow velocity 
can be neglected [S]. Thus the fluctuation in Teg may be 
a good representation of the fluctuation in temperature, 
(Z\Tkin = ‘^Teff). Being a thermodynamic quantity, the 
heat capacity is only connected to Tkin, and eqn.[2]can 
be rewritten as: 


1 {AT^.J (Z\T,^ff) 

7^2 ~ 7^2 

^ -^kin -^kin 


Event-by-event global temperature fluctuation has 
contributions from both the statistical and dynamical 
components. The statistical component arises due to 
limited number of particle production in heavy-ion col¬ 
lisions, and does not contribute to the measurement of 


heat capacity of the system. Statistical component of 
the fluctuation has been estimated by constructing ran¬ 
domly generated “synthetic” events, keeping the mul¬ 
tiplicity and pt distributions similar to those of the 
AMPT events. The (tox) and Teg distributions of the 
synthetic events are also shown in Fig. [S] The dynam¬ 
ical component of the global temperature fluctuation, 
obtained by subtracting the statistical component of 
the fluctuation from the total, has been estimated to be 
Z\Teg = 3.12 MeV. According to Ref. [31], Tkin = 0.095 
and (/3x) = 0.651 from combined Blast Wave hts. This 
corresponds to heat capacity of C = 927.1. 

Specific heat of the system is defined as heat ca¬ 
pacity per volume at the thermodynamic limit. For 
ideal gas, specific heat is (3/2)N, where N is the Avo- 
gadro’s number. For finite number of particles, the spe- 
cihc heat can be expressed as heat capacity per par¬ 
ticle. We define the specific heat (c„) as the heat ca¬ 
pacity per pion multiplicity within the central rapidity 
(—0.5 < y < 0.5) and 0 < (?) < tt). The c-u is found to 
be 0.633. Considering the variation of /3 t within 0.63 
and 0.67 [34], the value of Cy turns out to be within the 
range of 0.594 to 0.687. This is the Hrst estimation of 
Cy at LHC energies from model calculations. 


4 Local Temperature Fluctuations 

Local temperature fluctuations, which provide the amount 
of non-uniformity within a single event, are studied by 
dividing the available phase space into several y-(j) bins, 
and estimating the bin temperature (Tgin)- It helps to 
hnd local hotspots created during the initial energy den¬ 
sity, whether those have survived or died out. The bin 
temperatures are obtained using the similar prescrip¬ 
tion as above. The (tox) of pions are calculated within 
the y-cj) bin and toq < tox <1.5 GeV, and Tgin is eval¬ 
uated by using eqn. ®. The number of y-cj) bins has 
been chosen by taking into account the number of pi¬ 
ons in each bin so that fluctuations in the number of 
pions do not affect the estimation of tox- The amount 
of local fluctuation may vary depending on the number 
of y bins, which needs to be evaluated. 

For a given event, local temperature fluctuation in 
a given y-(f> bin is expressed as: 

Tbin = (Tbin — Teg)/Teg. (5) 

For each event, a fluctuation map in y-cj) phase space 
is constructed by plotting the corresponding values of 
Tbin- Fig.[4]shows the the temperature fluctuation map 
for a typical event in 6x6 bins in y-(j), where the fluctu¬ 
ations are represented by different colour pallets. This 
map gives a quantitative view of the local temperature 
fluctuation in the available phase space for an event. 
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Fig. 4 (Color online). Temperature fluctuation map in y- 
<f> bins for central Pb-Pb collisions at ^^snn = 2.76 TeV using 
the AMPT model. For each y-cj) bin, fluctuation is expressed 
as (Tbin — Teff)/Teff, the deviation of the bin temperature 
to the event temperature. The colour palettes indicate the 
magnitude of fluctuations. 


The map shows several hot (red) as well as cold (blue) 
zones, and zones with average (green) fluctuation through¬ 
out the phase space. It is to be seen whether the hot and 
cold zones have their origin from the extreme regions of 
phase space that existed during the early stages of the 
reaction. 

For a single event, the amount of local tempera¬ 
ture fluctuation is quantified by the ratio of RMS to 
the mean of the Tbin distribution. Local fluctuations 
in (mx) and corresponding temperature for each event 
have been evaluated and their event-by-event distribu¬ 
tions are plotted in Fig. [SJ The left and right panels of 
the figure show the distributions from (tot) and Tbin, 
respectively. Mean value of the local temperature fluc¬ 
tuation for the event sample is 12.98% for 6x6 bins in 
central rapidity. Statistical component of the local fluc¬ 
tuations has been extracted from the “synthetic” events 
as discussed earlier. The mean value corresponding to 
the statistical component of the local temperature fluc¬ 
tuation has been estimated to be 10.79% as shown in 
Fig. [S] The average dynamical local fluctuation, ex¬ 
tracted after subtracting the statistical component, for 
AMPT events is 7.2%. The non-zero value of local tem¬ 
perature fluctuation may imply that these might have 
the contributions from the early state fluctuations. 


5 Remarks 

Extraction of temperature fluctuations from experimen¬ 
tal data maybe affected by some of the effects. Event 
plane orientation is one such effects, which needs to 
be taken care of, especially for non-central events. Eor 



AmT/niT 


Fig. 5 (Color online). Event-by-event local temperature fluc¬ 
tuations, obtained from 6x6 y-cj) bins in central rapidity and 
full azimuth for central Pb-Pb collisions at .^/snn = 2.76 TeV 
using the AMPT event generator. 


the present study using AMPT, the events are event 
plane oriented. Fourier decomposition of the momen¬ 
tum distributions in the transverse plane yields a (j)- 
independent, axially symmetric radial flow component 
and a 0-dependent part containing the anisotropic flow 
coefficients. For most central collisions, radial flow re¬ 
mains similar for all the events and the anisotropic flow 
components do not affect the slope of the mx distri¬ 
bution. Final state effects, such as resonance decay, 
and hadronic rescattering tend to make the tox spec¬ 
tra softer, and so choice of the tox window has to 
be made in order to minimize such effects. Although 
present analysis uses charged pions, species dependence 
of temperature fluctuations may provide extra infor¬ 
mation regarding their freeze-out hyper surfaces as the 
particle production mechanisms of mesons, baryons and 
strange particles are different. This study may shed 
light on whether the origin of the temperature fluc¬ 
tuations are solely due to initial state fluctuations or 
includes final state effect. Viscosity tends to dilute the 
fluctuations. The SM version of AMPT includes the ef¬ 
fect of viscosity {r]/s ^ 0.15 at T=436 MeV [3M38] ). 
Further analysis using a viscous hydrodynamic models 
will be more realistic for this study. 

Local temperature fluctuation map for each event, 
as shown in Fig. |4j has a striking similarity to the fluc¬ 
tuation map of the cosmic microwave background radi¬ 
ation (CMBR) [31]. Fluctuation analysis of CMBR fluc¬ 
tuation map confirms the Big Bang evolution, inflation 
and provides information regarding the early Universe. 
The study of higher order moments using the maps may 
give access to various thermodynamic parameters at the 
early stages of the evolving system. Similarly, the fluc¬ 
tuation maps of heavy-ion collisions may form the basis 
of power spectrum analysis [4011311142]. Access to large 
number of events in heavy-ion collisions compared to 
single event analysis in CMBR may have definite ad¬ 
vantage which can be utilized to our advantage in or- 
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der to gain access to conditions that prevailed at the 
primordial state. 


6 Summary 

Event-by-event temperature fluctuations over full phase 
space as well as local phase space bins have been pro¬ 
posed to characterize the hot and dense system pro¬ 
duced in heavy-ion collisions at ultra-relativistic ener¬ 
gies. The global temperature fluctuations provide the 
heat capacity as well as specific heat of the system, 
whereas the observation of local fluctuations would im¬ 
ply the presence of fluctuations at early stages of the 
collision. Relativistic hydrodynamic calculations have 
been used to understand the evolution of e and T fluc¬ 
tuations. It shows that the system exhibits fiercely large 
fluctuations at early times, which diminishes with the 
elapse of time. 

The feasibility of studying temperature fluctuations 
in Pb-Pb collisions at ^snn= 2.76 TeV has been demon¬ 
strated by using simulated events from the AMPT model. 
Temperatures are extracted from (mx) of charged pi¬ 
ous. The global fluctuation in the event temperature 
has been extracted and used to calculate the heat ca¬ 
pacity as C = 927.1. The specific heat per pion has 
been extracted to be c„ = 0.633 ± 0.054 at the LHC 
energy, which is consistent with the results at RHIC 
energies m- At the LHC, the phase transition is ex¬ 
pected to be a cross over. Thus the transition may not 
take place at a unique point in the phase diagram due 
to the evolution of entropy fluctuation created at the 
initial state, which reflects into energy density and tem¬ 
perature fluctuations. Thus, the estimation of the heat 
capacity helps in understanding the nature of the phase 
transition. The variation of Cy as a function of center-of- 
mass beam energy for Au-|-Au collisions at RHIC may 
provide an effective tool for locating the QCD critical 
point [^131) . At LHC energies, it is possible to extract 
local temperatures over small phase space bins in cen¬ 
tral rapidity. Extraction of Local temperature fluctua¬ 
tion complements the global event temperature fluctu¬ 
ation as the origin of local fluctuations may be of pri¬ 
mordial in nature. For the AMPT, local temperature 
fluctuations over small phase bins within each event 
have been extracted. The amount of local fluctuation is 
7.2% for 6x6 bins in central rapidity. The observation 
of the non-zero local fluctuation may imply that a part 
of the initial fluctuations might have survived till freeze- 
out. The present study of global and local temperature 
fluctuations in conjunction with theoretical model cal¬ 
culations open up new avenues for characterizing the 
heavy-ion collisions. 
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